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ABSTRACT 


Seismic  data,  was  digitally  collected  from  the  recently  extended 
short-period  array  of  the  Tonto  Forest  Observatory  during  the  summer  of  1968 
and  the  winter  of  1969.  An  analysis  of  the  data  was  made  to  characterize  the 
ambient  noise  field  from  each  period  and  to  determine  time-and  space- station  - 
arity  of  the  fields.  The  effectiveness  of  multi-channel  filter  processing  on 
the  extended  array  was  also  measured. 

Several  contributors  to  the  ambient  noise  field  could  be  defined 
in  both  of  the  periods  analyzed  and  appeared  to  be  both  time^lind  space -station¬ 
ary.  Other  contributors,  due  primarily  to  atmospheric  storm  or  pressure 
activity,  were  not  stationary  and  could  be  related  to  its  generating  source. 

Multi-channel  filtering  of  the  extended  array  could  not  provide 
any  significant  improvement  over  simple  beamsteering  at  frequencies  of  interest 
(greater  than  .8  Hz).  At  lower  frequencies,  the  amount  of  improvement  obtained 
by  MCF  processing  relative  to  beam  steering  was  directly  related  to  the  increase 
in  noise  level  due  to  the  non -stationary  contributor  of  the  noise  field. 
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SECTION  I 


INTRODUCTION 

Two  data  collection  efforts  were  conducted  at  TFO  for  the  pur¬ 
pose  of  recording  high  quality  digital  data  from  the  recently  extended  short- 
period  array  and  the  three  component  long-period  array.  The  data  collecfed 
were  classified  as  'summer'  data,  consisting  of  the  data  collected  during  the 
pe-iod  of  6  June  to  21  June  1968,  and  as  'winter'  data,  which  was  recorded 
during  the  period  of  27  January  to  4  March  19o9. 

This  report  covers  the  TFO  extended  short-period  array  noise 
study.  The  primary  objectives  of  the  task  were  to  compare  the  'summer'  and 
'winter'  noise  field  characteristics,  to  estimate  the  time  and  spacial  stationarity 
of  these  noise  fields,  and  to  establish  the  effectiveness  of  multichannel  process¬ 
ing  using  the  extended  array. 

The  analysis  of  the  ambient  noise  field  and  the  effectiveness  of 
multichannel  processing  on  the  array  outputs  was  conducted  almost  exclusively 
in  the  frequency  domain  using  crosspower  matrices  developed  from  noise 
samples  selected  from  each  of  the  recording  periods.  A  summary  of  the  analysis 
and  conclusion  are  presented  in  the  next  section.  A  detailed  description  of  the 
processing  and  analysis  of  the  data  is  presented  in  following  sections. 
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v  SECTION  II 

» 

•  *  *  SUMMARY  AND  CONCLUSIONS 

The  ambient  noise  field  of  the  extended  short-period  array  at 
TFO  was  characterized  from  a  spectral  analysis  of  data  samples  selected 
from  each  of  the  recording  periods.  The  time-  and  space -stationary  of  the 
coherent  noise  energy  was  measured  from  absolufjj.-'Single  channel  power- 


density  spectra  and  from  high  resolution  wavenumber  spectra.  The  effective- 

$ 

ness  of  multichannel  processing  on  the  extended  array  was  determined  from 

| 

the  design  and  Evaluation  of  MCFs  on  crosspower  matrices  of  data  from  each 

'  A 

of  the  periods  being  analyzed. 

A.  ANALYSIS  OF  THE  AMBIENT  NOISE  FIELD 


Absolute  level  of  the  ambient  noise  field  was  determined  from 
single -channel  power  density  spo.  a,  computed  by  the  Markov  technique. 
Significant  changes  in  the  absolute  level  of  the  power  frequency  energy  were 

evident.  The  'winter1  samples  showed  a  general  increase  in  the  absolute  power 

'  1  \ 

level  especially  at  frequencies  below  .25  Hz.  The  increases  in  the  absolute 
levels  of  the  power  spectra  appeared  to  be  directly  related  to  the  amount  and 

intensity,  of  atmospheric  pressure  activity  observed  from  weather  maps.  A 

\ 

significant  increase  in  energy  is  observed  in  the  area  of  .  12  Hz  during  a 
'winter'  period  of  very  intense  atmospheric  disturbance  off  the  Alaskan  Coast. 


The  spacial  organization  of  the  coherent  noise  energy  during 
the  winter  and  summer  period  was  determined  from  high  resolution  frequency- 
waveftumber  spectra  computed  for  each  of  the  periods. 


Surface  mode  energy  from  the  Northeast  is  observed  during  both 
periods.  This  energy  is  observed  with  the  fundamental  surface  node  velocity 
of  approximately  3.2  km/second  at  frequencies  below  .25  Hz.  At  higher  fre¬ 
quencies  this  energy  appears  to  travel  at  a  higher  velocity  of  approximately  3.6 
km/second  indicating  the  existence  of  higher-order  surface  mode  energy. 

Burg,  J.  P.  ,  "Maximum  Entropy  Spectral  Analysis"  Society  of  Exploration 
Geophysicists _ _ 
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At  frequencies  of  .4  Hz  and  greater  the  spacial  organization  of 
the  coherent  surface  mode  energy  becomes  isotropic  in  nature.  At  these  fre¬ 
quencies  a  significant  amount  of  P-wave  energy  from  the  South,  with  very  high 
apparent  horizontal  velocities  (40  km/second  or  greater)  are  consistently 
observed.  Secondary  P-wave  contributions  from  the  Northeast  and  Northwest 
with  slightly  lower  apparent  horizontal  velocities  are  also  observed. 

The  most  significant  difference  in  the  spacial  organization  of  the 
summer  and  winter  ambient  noise  is  the  addition  of  noise  from  the  Northwest 
as  observed  in  the  winter  noise  sample.  At  .  12  Hz,  the  surface  mode  energy 
from  this  direction  dominates  the  spectra.  However,  as  the  frequencies  increase, 
the  surface  mode  energy  contribution  from  this  direction  rapidly  decreases.  The 
P-wave  energy  from  the  Northwest  continues  to  be  a  significant  contributor  to  the 
noise  field  up  to  frequencies  of  approximately  .  5  Hz.  Analysis  of  the  measured 
apparent  horizontal  velocity  of  the  P-wave  energy  from  the  Northwest  indicates 
the  Alaskan  storm  to  be  the  generating  source. 

A  slight  increase  in  the  surface  mode  energy  from  the  South  and 
Southwest  is  also  observed  in  the  'winter'  data.  In  general  the  increase  in  the 
surface  mode  energy  relates  very  well  to  the  pressure  activity  in  the  correspond¬ 
ing  direction. 

High-resolution  K-line  spectra  were  computed  to  determine  the 
relative  contribution  of  each  noise  source  to  the  total  ambient  noise  power.  The 
percent  contributed  as  P-wave,  directional  surface  mode,  and  isotropic  surface 
mode  is  illustrated  in  Figure  IV- 18. 

Based  on  the  spectral  analysis  of  the  noise  from  the  two  time 
periods, the  spatial  organization  of  certain  components  of  the  coherent  noise 
field,  surface-mode  energy  from  the  Northeast  and  P-wave  from  the  South, 
are  stationary.  However,  the  other  surface-mode  and  P-wave  components, 
which  are  highly  influenced  by  seasonal  and  daily  atmospheric  variations  or 
distrubances ,  are  not  time  or  space  stationary. 
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B.  ANALYSIS  OF  MCF  PER  FOP  MANCE 


A  set  of  multichannel  filters  were  designed  from  a  crosspower 

/ 

matrix  computed  from  a  long  noise  sample  fron,  each  of  the  time  periods.  The 
signal-to-noise  in^provement  was  determined  from  a  frequency  domain  appli¬ 
cation  of  the  MCF  to  the  noise  and  signal  crosspower  matrices  used  in  the 
design.  Results  were  compared  to  a  similar  applicaton  of  a  simple  beamsteer 
(straight  summation)  process  to  the  same  matrices. 

Using  34  elements  of  the  extend  short-period  array  a  S/N  improve 
ment  of  23  db  was  obtained  by  MCF  processing  at  the  frequency  of  .  12  Hz  on  the 
'winter*  data.  Wavenumber  analysis  of  MCF  response  shows  this  to  be  due 
primarily  to  the  rejection  of  the  surface-mode  energy  fro  r»  the  Northwest. 

In  general,  some  S/N  improvement  can  be  obtained  by  MCF  pro¬ 
cessing  over  that  of  beamsteering  at  frequencies  below  .  75  Hz.  On  the  'summer* 
data,  S/N  improvements  of  3  to  6  db  alpove  that  of  the  beamsteering  processing 
was  obtained  by  MCF  processing.  Additional  gain  is  obtained  on  noise  samples 
which  have  an  increase  in  coherent  energy  at  these  frequencies  due  to  the  pre¬ 
viously  discussed  atmospheric  activity.  However,  since  the  generating  source  is 
not  considered  as  space  stationary,  on  line  processing  with  fixed  MCFs  is  not 
expected  to  provide  any  significant  improvement.  At  frequencies  above  .8  Hz, 
MCF  processing  of  the  extended  array  would  not  provide  any  significant  signal- 
to-noise  improvement  over  beamsteering  even  under  stationary  conditions. 
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SECTION  HI 

PREPARATION  OF  DATA  FOR  ANALYSIS 

Large  quantities  of  seismic  data  were  collected  from  the  extended 
short-period  array  on  a  Digital  Field  System  (DFS)  operated  by  Texas  Instru¬ 
ments'  engineers  during  two  recording  efforts.  The  first  recording  session 
covered  the  period  of  6  June  1968  to  21  June  1968.  Data  from  this  period  will 
be  designed  as  'summer'  data.  The  second  recording  period,  during  which 
'winter'  data  were  recorded  covered  the  period  from  27  January  1969  to  3  March 
1969.  This  section  outlines  the  pre-processing  performed  on  noise  samples 
selected  from  the  recorded  data  for  extensive  analysis.  For  reference.  Figure 
HI- 1  shows  the  relative  locations  of  the  extended  short-period  array  elements 
from  which  the  data  was  collected. 

A.  EDITING  OF  DATA 

The  'summer'  data  was  recorded  during  a  period  when  severe 
spiking  problems  on  array  elements  were  being  experienced  on  the  recently 
extended  array.  Several  segments  of  the  data,  covering  approximately  two 
hours,  were  selected  and  edited  from  the  field  tapes.  One  of  these  samples 
recorded  on  19  June  1968,  containing  33  channels  of  seismic  data  (elements 
Z32,  Z6,  Z18,  and  Z20  were  inoperative)  was  then  processed  with  a  program 
to  detect  and  remove  the  undesired  spikes  from  the  data. 

The  data  was  recorded  by  the  DFS  with  a  sampling  interval  of 
.096  seconds.  Since  the  noise  energy  at  the  higher  frequencies  (greater  than 
1.5  Hz)  was  expected  to  be  random  across  the  extended  array  due  to  the  large 
inte r element  spacing,  the  data  was  high-cut  ianti-aliased)  filtered  and  resampled 
to  provide  a  sampling  interval  of  .  288  seconds  and  a  frequency  band  of  0  -  1.  736 
Hz.  The  data  was  also  pre-whitened  as  shown  by  Figure  IH-2. 
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Figure  III-2.  Spectra  Illustrating  Effects  of  Anti-Alias  and 
Pre-Whitening  Filter. 

Similarly,  several  segments  of  data  were  edited  from  the  'winter' 
data  tapes.  One  of  the  noise  samples  recorded  on  27  february  1969  during  which 
intense  low  pressure  conditions  were  observed  from  weather  maps,  was  selected 
for  a  detailed  analysis.  However,  due  to  bad  spots  in  the  data  such  as  missing 
elements  and  severe  spiking, only  50  minutes  of  data  from  any  recording  was 
acceptable  for  analysis.  All  but  three  (Z19,  Z22,  and  Z31)  of  the  37  array 
elements  were  operating  satisfactorily.  The  data  was  anti-aliased  filtered  and 
resampled  to  provide  the  .  288  second  sampling  interval. 

B.  DEVELOPMENT  OF  X-POWER  MATRICES 

The  fast  Cooley- Turkey  transform  method  was  used  to  obtain  the 
frequency  domain  representation  of  the  'summer'  and  'winter'  noise  samples. 

The  crosspower  matrices  were  developed  for  each  of  the  noise  samples  as 
follows. 

HI- 3  service*  group 
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1.  'Summer'  Data 


The  time  traces  of  the  'summer'  data  was  transformed  using 
gates  of  512  points  of  data  at  a  time.  A  total  of  45  gates  were  transformed. 
Crosspower  matrices  for  the  33  channels  of  data  were  computed  from  the  trans¬ 
forms  for  each  of  the  257  frequencies  covering  the  band  0.0  to  1.736  Hz  in 
.00678  Hz  increments.  The  matrices  were  then  stacked  for  all  the  transform 
gates  and  smoothed  across  frequencies  by  an  equal  weighting  function  for  five 
adjacent  frequencies.  Every  fifth  frequency  was  saved  for  further  analysis  re¬ 
sulting  in  a  set  of  51  matrices  for  the  0-1.736  Hz  band  at  a  frequency  incre¬ 
ment  of  .0339  Hz  between  matrices.  The  spectral  window  due  to  the  transform 
gate  and  smoothing  described  is  shown  in  Figure  IH-3. 

2.  'Winter'  Data 

The  'winter'  time  traces  were  transformed  in  256  point  segments 
resulting  in  129  complex  frequency  values  for  the  0  to  1.736  Hz  frequency  band. 
A  total  of  40  segments  were  transformed  for  the  34  channels  of  data.  Cross¬ 
power  matrices  were  computed  from  the  transforms  and  stacked.  In  this  case, 
smoothing  of  the  matrices  was  performed  across  three  adjacent  frequencies 
with  equal  weighting  of  matrices.  This  resulted  in  a  set  of  43  matrices  covering 
the  frequency  band  of  the  data  in  .041333  Hz  increments.  The  effective  spectral 
window  for  this  data  is  also  shown  in  Figure  III- 3. 
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Figure  III-3.  Spectra  Window  of  the  Crosspower  Matrix  Computed  from  the  TFO  Noise  Samples 


SECTION  IV 


ANA T  YSIS  OF  AMBIENT  NOISE  FIELDS 

The  characteristics  of  the  'winter'  and  'summer'  noise  field 
were  determined  and  compared  through  an  analysis  of  single  channel-power 
density  spectra,  frequency-wavenumber  spectra  and  K-line  spectra  for  the 
extended  short-period  array.  The  time-  ai  d  space- stationarity  of  the  noise 
field  were  estimated  from  the  noise  samples  analyzed. 

A.  SINGLE -CHANNEL  POWER -DENSITY  ANALYSIS 

Power  density  spectra  were  estimated,  using  the  Markov  Spectral 
Estimation  Technique  for  several  channels  from  each  of  the  noise  samples 
studied. 

The  daily  1.0  Hz  calibration  of  the  array  sensors  were  recorded 
by  the  DFS  and  used  in  this  study  to  establish  absolute  levels  for  the  spectra. 

The  resulting  spectra  for  the  'winter'  and  'summer'  data  are  shown  in  Figures 
IV-1  and  IV- 2  respectively. 

A  comparison  of  the  spectra  show  an  increase  of  approximately 
6  db  in  noise  power  for  the  'winter'  sample  across  all  frequencies  shown.  A 
larger  power  increase  of  approximately  12  db  is  observed  in  the  area  of  .25  Hz. 

During  the  period  which  the  'winter'  data  was  recorded,  intense 
lows  occurred  at  several  locations  along  the  coast  lines.  The  strongest  of  these 
were  situated  off  the  west  coast  of  Alaska.  Others  existed  along  the  west  coast 
of  continental  U.S.  ,  west  coast  of  Mexico,  Gulf  of  Mexico,  in  the  Baffin  Bay, 
east  coast  of  Greenland,  and  along  the  eastern  coast  line  of  U.S.  near  the  New 
England  states.  The  'summer'  sample  used  in  the  analysis  was  recorded  during 
a  period  of  relative  little  storm  activity.  However,  another  'summer'  noise 
sample  was  selected  for  the  single-channeLspectral  analysis,  which  was  recorded 
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Figure  IV- 1.  Single-Channel  Power  Density  Spectra  From  'winter*  Record¬ 
ing  Period 


IV -2 


services  group 


TO  10  my'At  AT  1  OKI  TO  l  0  mu '  AT  l  i|  *  ?0  Li  *^*40  AT  10*  TO  l.t  AT  1 1  * 


1 


serves  group 


i 


during  a  period  when  some  intense  atmospheric  activity  was  occurring  off  the 
Newfoundland  coast.  This  sample  10  June,  will  be  classified  as  'noisy  summer1 
data.  A  comparison  of  the  spectra  from  the  'winter'  and  two  'summer'  data  sam^ 
pies  is  presented  in  Figure  IV-3.  The  'noisy  summer'  sample,  has  a  compara¬ 
ble  noise  level  to  the  'winter'  data  in  the  .  3  to  1.  0  Hz  frequency  band.  It  does 
not  show  the  large  increase  in  noise  energy  below  .  3  Hz  as  observed  in  the 
'winter'  data.  The  following  wavenumber  analysis  will  show  this  increase  in 
energy  below  .  3  Hz  in  the  'winter'  noise  to  be  due  primarily  to  the  intense  low 
pressure  activity  off  the  Alaskan  coast. 

B.  HIGH-RESOLUTION  FREQUENCY  WAVENUMBER  ANALYSIS 

it- 

The  crosspower  matrices,  <t>  ,  developed  in  Section  III  for  each 
of  the  noise  samples  were  used  to  compute  high -resolution  wavenumber  spectra 
at  frequencies  in  the  .  12  to  .6  Hz  band.  Spectra  from  higher  frequencies  become 
difficult  to  interprete  due  to  the  spacial  alaising  of  the  array.  The  method  used 
to  compute  the  spectra  is  described  by: 


_  2 

P(k)  = 
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NC  =  number  of  array  sensors, 
V  is  the  column  vector 
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-Comparison  of  'winter'  and  'summer'  -Single-/Channel  Power  Density  Spectra 


and  T  is  the  reference  vector  for  the  sene^r,  i.e.  all  components  are  zero 
i 

except  the  one  corresponding  to  the  sensor  being  referenced. 


where  a  is  a  positive  non-zero  value. 

•  v 

1.  'Slimmer'  Surface-Mode  Spectra 

High -resolution  wavenumber  spectra  computed  at  the  frequencies 
of.  169,  .203,  .236,  .270,  .304,  .338,  .372,  and  .  405  Hz  from  the  'summer* 
noise  matrices  are  shown  in  Figure  IV-4  through  IV-7.  The  edge  velocity  of 
each  of  the  spectra  is  2.  0  km/second.  A  dashed  circle  is  shown  on  the  spectra 
to  illustrate  the  approximate  velocity  of  surface  mode  energy. 

•  In  the  first  figure,  the  surface  mode  energy  from  the  Northwest 

•  1  ,.V 

appears  to  be  a  principal  contributor  to  the  noise  field  along  with  the  P-wave 
energy  from  the  South.  Measurements  of  the  surface  mode  peak  indicates  an 
apparent  velocity  across  the  array  of  3.  17  km/second,  A  secondary  surface 
mode  peak  at  approximately  the  same  velocity  is  obs*«!>ed  from  the  Southwest. 

As  the  frequencies  are  increased  the  surface  mode  peaks  rapidly 
’  "*  •  ¥  .  ' 

decrease  in  relative  power  and  the  p-wave  from  the  south  becomes  predominant. 

The  surface  mode  energy  becomes  very  isotropic  in  nature  with  a  measured 

velocity  of  3.  2  krn/second.  There  is  continued  evidence  of  directional  surface 

mode  energy  from  the  Northeast  at  the  higher  frequencies,  however,  the  apparent 

velocity  is  somewhat  higher  than  that  previously  measured  for  the  directional 

'contribution  and  that  of  the  isotropic  noise  indicating  that  it  is  traveling  as 

some  higher  order  mode  energy. 
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Figure  IV-5.  High  Resolution  Wavenumber  Spectra  of  .23  and  .27  Hz 
'summer'  Noise 
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The  energy  arriving  from  the  Northeast  is  believed  to  be  gen¬ 
erated  in  the  area  of  Newfoundland  and  the  eastern  coast  of  Canada  and  the  US 
based  on  the  spacial  organization  of  P-wave  energy  and  the  correlation  of  atmo¬ 
spheric  activity  in  this  region  to  increased  surface  mode  noise  from  this  direction. 
The  path  of  such  surface  mode  noise  would  transverse  the  Rocky  Mountain  region 
which  could  possibly  cause  the  higher  order  mode  energy  observed  at  the  higher 
frequencies.  That  is  the  higher-mode  excitation  might  have  been  effected  by  the 
fundamental  mode  impinging  on  the  discontinuity  of  the  mountainous  range. 

2.  'Summer'  P-Wave  Spectra 

At  the  higher  frequencies  the  resolution  of  the  array  becomes 
sufficient  to  show  the  P-wave  noise  field  consists  of  several  contributors.  In 
an  effort  to  analyze  the  P-wave  noise  in  detail,  wavenumber  spectra  with  an 
edge  velocity  of  10.  0  km/second  were  computed  and  are  shown  in  Figures  IV-8 
and  IV-9  for  the  frequencies  of  .  507,  .  541,  .  574,  and  .  608  Hz. 


As  observed  in  these  spectra  the  principle  p-wave  contribution 
is  from  the  south  with  an  apparent  horizontal  velocity  of  30  to  40  km/second. 
As  the  frequencies  increase,  the  direction  of  this  energy  shifts  slightly  from 
east  to  west. 


A  contributor  to  the  p-wave  field  is  also  shown  in  the  same  direc¬ 
tion  as  the  earlier  observed  surface  mode  energy  from  the  Northeast.  Weather 
maps  of  this  period  show  low  pressure  activity  near  the  southern  tip  of  Green¬ 
land.  Events  occurring  in  this  area  would  arrive  at  TFO  with  an  azimuth  of  31 
degrees  and  an  apparent  horizontal  velocity  of  15.  5  km/second.  The  average 
velocity  and  azimuth  of  the  Northeast  peak  in  Figure  IV-9  is  15  km/second  and 
40  degrees,  respectively,  indicating  gener  ating  source  to  be  in  general  area  of 
the  atmospheric  activity. 

Similarly,  the  p-wave  contribution  from  the  Northwest  can  be 
related  to  low  pressure  activity  below  Alaska.  The  azimuth  of  the  energy  peak 
in  the  wavenumber  spectra  is  measured  to  be  approximately  315  degrees  and 
the  velocity  range  is  11.  3  to  15.  8  km/second.  An  event  occurring  off  the 
Southern  coast  of  Alaska  in  the  area  of  the  storm  activity  would  arrive  at  TFO 
with  an  azimuth  of  323  degrees  and  a  velocity  of  12.4  km/second. 
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3. 


'Winter'  Surface-Mode  Spectra 


Analysis  of  the  single-channel  power  density  spectra  showed 
significant  energy  increase  at  very  low  frequencies  in  the  'winter'  sample. 

A  high  resolution  wavenumber  spectra  computed  at  the  frequency  of  .  12  Hz 
(Figure  IV-10)  indicates  this  increase  resulted  primarily  from  a  strong 
surface-mode  energy  contribution  from  the  Northwest  in  the  direction  of  the 
previously  mentioned  intense  low  pressure  disturbance  off  the  Western  coast 
of  Alaska.  However  the  high  resolution  technique  also  indicates  a  second 
peak  in  the  spectra  at  this  frequency  is  1  80°  to  the  previous  one  and  is  at  a 
somewhat  lower  velocity.  An  analysis  of  wavenumber  spectra  computed  by 
conventional  methods  were  compared  to  the  high  resolution  spectra  to  estab¬ 
lish  authenticity  for  each  of  the  observed  energy  lobes.  Conventional  spectra 
computed  at  .  1  2  Hz  gave  no  indication  of  noise  energy  from  the  direction  of 
this  second  lobe.  It  is  more  likely  due  to  effects  resulting  from  the  strong 
highly  coherent  noise  source  from  the  Northwest.  Similar  effects  have  been 
observed  in  wavenumber  responses  of  MCF's  designed  using  noise  statistics 
consisting  of  a  single  highly  coherent  noise  component.  At  the  higher  fre¬ 
quencies,  the  percentage  of  the  total  coherent  noise  energy  contributed  by  a 
single  source  is  significantly  less  than  that  observed  for  the  Northwest  com¬ 
ponent  at  .  I  2  Hz,  thereby  restricting  such  effects. 


IV- 14 


services  group 


Figure  IV- 10.  High  Resolution  Wavenumber  Spectra  of  the  .  12  Hz 
'winter'  Noise 


In  Figure  IV- 11  the  wavenumber  spectra  for  the  frequencies  of 
.  163  and  .203  Hz  are  shown,  in  which  additional  surface  mode  }obes,  similar 
to  those  observed  in  the  'summer'  sample,  are  evident.  The  multiple  surface 
mode  lobes  are  in  close  agreement  with  expectations  based  on  the  numerous 
atmospheric  disturbance  observed  from  weather  maps. 

At  the  higher  frequencies  of  .24  to  .41  Hz  (Figures  IV-12  and 
IV- 13),  the  generating  source  to  the  Northeast  of  the  array  becomes  the  princi¬ 
ple  contributor  to  the  surface-mode  energy  a$  previously  observed  in  the 
'summer'  data.  However,  the  principle  p-wave  contributor  for  the  'winter' 
noise  is  from  the  Northwest  ^Alaskan  activity)  contrary  to  the  predominant 
Southern  P-wave  contributor  in  the  'summer'  data. 
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4.  'Winter1  P-Wave  Spectra 

H:gh  resolution  wavenumbe."  spectra  were  computed  with  an  edge 
velocity  of  10  km/second  at  the  higher  frequencies  of  .41,  .49,  .57,  and  .61  Hz 

J 

to  determine  the  spacial  structure  of  the  'winter'  p-wave.  These  spectra  are 
shown  in  Figure  IV- 14  and  IV- 15. 

At  .41  Hz  the  predominant  contributor  of  p-wave  is  observed  to 
be  from  the  Northwest  and  a  secondary  contribution  from  the  South.  As  the 
frequencies  are  increased  the  Southern  contribution  continually  increases  in 
relative  energy  contribute  u  over  that  of  the  Northwest  source. 

The  average  velocity  and  azimuth  of  the  p-wave  energy  from  the 
Northwest  was  measured  from  the  four  spectra  to  be  15.  3  km/second  and  320 
degrees.  P-wave  energy  originating1  in  the  area  of  the  atmospheric  activity 
near  the  Western  coast  of  Alaska  would  arrive  at  TFO  with  an  apparent  hori¬ 
zontal  velocity  of  13.  9  km/second  and  an  azimuth  of^  326  degrees. 

As  in  the  'summer'  data,  the  p-wave  energy  due  to  the  Southern 

contribution  has  a  very  high  apparent  horizontal  velocity  of  40  km/second  or 

higher.  A  possible  generating  source  of  the  p-wave  energy  from  the  South  could 

.  * 

not  be  defined. 
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C.  K-LINE  SPECTRAL  ANALYSIS 


The  1 -dimensional  wavenumber  (K-line)  spectra  presented  and 

analyzed  in  this  section  can  be  interpreted  as  the  physical  projections  of  the 

2-dimensional  power  density  wavenumber  spectra,  similar  to  that  discussed 

in  previous  subsection,  onto  axes  which  are  parallel  to  some  line  of  array 

sensor.  The  3-lines  of  approximately  equally- spaced  sensors  used  in  this 

study  is  illustrated  as  arms  in  Figure  III-l.  The  resulting  K-line  spectra  for 

each  line  gives  the  power  density  of  the  ambient  noise  as  a  function  of  its 

apparent  wavenumber  along  that  line.  A  more  detailed  discussion  of  the  K-line 

spectral  technique  has  been  published  in  a  special  report  on  K-line  spectral 

$ 

analysis  using  the  TFO  cross-array. 


udoic  input  Uc 

ix,  £  N.j(f)  J  , 


^  The  basic  input  data  for  calculating  the  K-line  spectra  in  the 
cros3power  matrix,  £  N„(f)  J  ,  where  f  is  frequency  and  i  and  j  range  over  the 
sensors  along  the  line.  The  crosspower  matrix  developed  in  Sectidn  III  for 
the  entire  extended  array  was  used  to  obtain  such  a  matrix  for  each  line  of 
sensors.  The  matrix  from  all  parallel  lines  of  sensor  were  averaged  to  obtain 

a  single  matrix  for  each  of  three  directions  for  which  a  K-line  spectra  were 

«<  • 

computed.  The  average  matrix  was  then  normalized  so  that  the  main  diagonal 
elements  were  unity. 

The  approximate  spacing  of  the  sensors  forming  the  lines  was 
5  kilometers.  The  resulting  foldover  wavenumber  of  the  spectra  was  .  1 
cycles /km. 

The  K-line  spectra  computed  for  the  three  lines  from  all  the 
matrices  generated  from  the  'winter1  and  'summer'  data  in  the  fre,quency 
band  .  1  to  .  5  Hz  are  in  Figure  IV-16  and  IV- 17  respectively.  Each  plot  consists 


Array  Research,  Analysis  of  K-line  Waverifember  Spectra  from  the  TFO 
Long- Noise  Sample,  Special  Report  No.  23.  -> 


of  a  wavenumber  power  density  spectrum,  an  integrated  wavenumber  power 
density  function,  and  a  discrete  function  marked  by  X's  which  show  the  frac¬ 
tional  part  that  is  unpredictable  during  an  attempt  to  predict  the  output  of  the 
next  sensor  in  a  line  from  that  of  the  preceding  ones  of  the  line.  The  0-db 
level  of  the  right  hand  scale  for  the  power  density  spectrum  is  the  average 
value  of  the  spectrum.  The  fractional  scale  for  the  integrated  spectra  and  the 
prediction  error  is  on  the  right  hand  side.  Two  sets  of  vertical  dashed  lines 
are  shown  to  indicate  the  approximate  velocity  of  surface  mode  energy  either 
traveling  along  the  line  or  isotropic  in  nature  (outer  set)  and  to  show  the 

minimum  apparent  velocity,  of  p-wave  noise  (inner  set).  The  inner  set  of 

f> 

vertical  lines  are  not  to  be  confused  as  an  envelope  for  p-wave  energy  since 
directional  surface  mode  energy  could  also  fall  inside  them  depending  on  the 
angle  of  incidence  of  this  energy  to  the  line  of  sensors. 

A  simultaneous  analysis  of  the  three  K-line  spectra  at  each 
frequency  confirms  the  spacial  organization  of  the  ambient  noise  previously 
observed  in  the  2-dimensional  wavenumber  analysis.  The  complexity  and 
higher  prediction  of  the  'winter'  spectra  compared  to  the  'summer'  spectra  is 
due  to  the  multiple  noise  sources  in  the  'winter'  sample  as  previously 
described. 


*■ 
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PERCENT  OF  PREDICTION  ERROR 
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In  the  previous  analysis,  the  direction  of  several  definite 
contributors  to  the  ambient  noise  field  could  be  determined.  The  integraded 
wavenumber  power  density  function  was  used  in  this  study  to  estimate  more 
accurately  as  a  function  of  frequency  the  percentage  each  contribution  is  mak¬ 
ing  to  the  total  'summer'  and  'winter'  noise  level.  The  contributors  considered 
consist  of  4  directional  surface  mode  sources,  isotropic  surface  mode,  and 
P-wave  as  illustrated  in  left  most  part  of  Figure  IV-18.  The  P-wave  con¬ 
tribution  consists  of  any  coherent  energy  with  an  apparent  velocity  of  8  km/ sec 
or  greater  and  would  be  contained  inside  the  8  km/sec  circle  of  the  wave- 
number  space  represented  by  the  figure.  The  surface  mode  contributors  would 
be  represented  in  wavenumber  space  along  the  circle  representing  the  surface 
mode  velocity  of  3.2  km/sec.  The  percent  of  th-  total  noise  field  which  each 
contributed  in  the  two  noise  samples  analyzed  are  illustrated  on  the  right  of 
Figure  IV-18.  The  most  significant  difference  observed  between  the  two  noise 
samoles  is  the  13rge  percentage  of  surface  mode  energy  in  the  'winter'  noise 
at  frequencies  below  .  24  Hz.  This  appears  to  be  primarily  due  to  the  energy 
from  the  direction  of  the  Alaskan  storm  velocity  previously  discussed  which 
significantly  increases  at  lower  frequencies.  Such  a  generating  source  of 
noise  should  be  of  significant  concern  for  longer  period  analysis  since  the 
percent  of  contribution  is  rapidly  increasing  for  the  longer-period  energy. 

A  similar  observation  can  be  made  of  the  surface-mode  energy  from  the 
South. 

/  ' 
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Figure. IV- 18.  Distribution  of  Ambient  Noise  Power  for  Identified  Sou 
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MU  i_/T  I  CHANNEL  FILTER  DESIGN  AND  EVALUATION 


The  desirability  of  multichannel  filtering  for  the  extended  short 
period  array  was  determined  by  designing  and  evaluating  signal  extraction 
fillers  from  the  crosspower  matrices  developed  in  Section  III.  The  perform- 
ance  of  the  multichannel  filters  were  compared  to  that  c,t  the  simple  beam¬ 
steering  process.  ,  . 

A.  MCF  DESIGN 

A  set  of  optimum  Weiner  multichannel  filters  weredesigned  in 
the  frequency  domain  from  each  of  the  generated  cros  jpower  matrices  develop¬ 
ed  fqr  the  'winter'  and  'summer'  data.  The  set  of  MCF?  from  the  'winte 
data,- denoted  as  MCF  W,  consist  of  3  MCF  designs,  each  using  a  subset  of 
the  34  array  elements.  The  first  'Winter'  MCF  (MCF  W-7)  was  designed  on 
the  seven  center  elements  (center  seismometer  Zl)  plus  the  intermost  ring  of 
seismometers)  of  the  extended  array.  MGF  W- 18  was  designed  on  thd  same 

elements  as  MCF’  W-7  plus  11  elements  of  the  next  ring  of  sensors  (sensors 

*  . 

Z19  and  Z22  were  not  operating  properly).  The  third  MCF*  MCF  W-34,  ,  | 

e  _  ^  0 

includes  all  of  the  elements,  except  Z19,  Z22,  and  Z3T.  .  — , 

t  ’ 

c  o  \  c  0 

An  equivalent  set  of  MCFs  developed  from  the  'summer'  data 
is  denoted  as  MCF  S-6,  MCF  S-  17,  and  Mcf  S-33.  In  this  case  only  6,  17, 

i.  '  •  •  v  - 

and  33  elements  re  spectively  we^re  used  due  to  one  of  the  elements  (Z-6)  on 
the  inner  ring  being  bad.  Figure  III-  1  is  shown  to  ihdicate  the  elements 
actually  used  in  each  of  the  designs. 

The  signal  model  used  in  the  MCF  designs  represented  waVes 

•  ;  •  -  T-  — ;  1-  * 

vertically  incident  to  the  extended  array.  The  signal  amplitude  on  each 

'  .  ..  >  ••  •  . 

channel  was  assumed* to  have  a  constant  ratio  (1)  to  the  noise  on  that  channel 


at  each  frequency.  J.h'e  signal  vector  V,  is  thus  defined  such  V  .  =yN.. where 
•  y  a  ^  ' 

N  .  is  an  element  of  the  noise  crosspower  matrix,  N.‘  The  signal  matrix  is 

defined  as  S<-  V  V^.  ?  ‘  -  - 
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The  equation  for  the  Wiener  filter  responses,  Z,  to  enhapce 

k  ’  1 

the  signal  with  minimum  mean  square  error, 

4b 

[N  +  S]  Z  =  (s/n)  V  where  (s/n)  =  4.  0, 

i 

was  solved  at  each  of  the  desired  frequencies. 

B.  EVALUATION  OF  VCFs 

The  designed  MCFs  were  evaluated  in  the  frequency  domain 

thrqugh  direct  applications  of  the  filter  response,  2?»  to  the  noise  matrix 

\  , 

from  whicih  they  were  designed  and  through  an  analysis  of  wavenumber 
responses  at  selected  frequencies.  The  S/N  improvement  obtained  in  the 
application  was  compared  to  that  obtained  by  simple  beamsteering  on  the 
same  data.  '  0  '  *  *  ; 

1.  Signal- to- Noise  Improvements  J 

v  The  signal- to- noise  improvement  for  each  of  the  MCFs  was 
obtained  by  taking  ratio  of  the  output  signal- to- noise  ratio,  S  /N  ,  where 

\  s  =  |vTzl  2  *  and 


where 


Z  NZ  ,  to  the  input  signal-to-noise  ratio,  S./N., 

»  ’  ii 


‘  S./N.  =  1.  / 

ii  ••  / 

c 

«  The  signal-to-noise  improvement  for,  the  set  of  'winter1  and 

'summer'  MOFs  is  presented  in  Figure  V-l  for  the  frequency  band  of  0.0  to 

<3  0  ' 

1.25  Hzt  A  comparison  of  the  two  sets  show  that  the  MCF  W  filters  obtain  a 
significantly  higher  S/N  improvement  below  0*75  . Hz  than,  the  equivalent  MCF  S 
filters,  as  expedtjsd  due  to  the  previously  observed  increase  of  surface  mode 

-ii  v‘  ,  0 

1  <a 

energy  in  the  'winter'  data  sample.  A  maximum  of  1 1  db  more  S/N  improve- 
ment  was  obtained  by  MCF  W-34  than  by  MCF  S-33  at  the  frequency  of  0.  12  Hz. 

’Similarly,  S/N  improvements  were  obtained  for  a  straight 
summation  process  oh  each  of  the  subset  of  elements.  A  ratio  of  the  MCF 

*  ,  *  ^  -  _i  ’  J  - 

improvement  to  this  summation  process  is  shown.in  Figure  V-2.  On  the 

.  ■  V_2  3  '  service®  group 


%  "  i1' 


TFO  SUMMER  DATA 


sorvicos  group 


’summer1  data  less  than  2.  0  db  gain  is  obtained  by  MCF  processing  at 

frequencies  higher  than  0.  8  Hzf  A  maximum  of  5.  6  db  S/N  improvement  over 

a  summation  process  was  obtained  at  0.  12^Hz  by  the  MCF  S-33.  The  MCF  W 

set  of  filters  show  significantly  more  improvement  over  the  summation.  This 

is  due  primarily  to  the  surface  mode  energy,  however,  since  less  statistics 

were  used  in  the  ’winter'  sample  than  the  'summer'  data  some  of  the  gain  is 

* 

false  due  to  overdesign  of  the  MCF.  Based  on  a  study  by  John  Burg,  et  al, 
resu  ts  from  evaluating  the  34  channel  filter  MCF  W-34  on  the  statistics  from 
which  it  was  designed  will  be  biased  approximately  1.4  db  over  the  true  opti¬ 
mum  Wiener  filter  performance.  (122  samples,  34  channels).  For  the  18 
channel  case,  MCF  W-18  the  evaluation  will  be  biased  only  .75  db.  For  the 
'summer'  data  230  effective  segments  of  data  was  used  to  generate  the  matrices 
thereby  resulting  in  a  33  channel  MCF  with  a  bias  of  less  than  .  7  db  relative  to 
the  Wiener  with  the  true  filter. 


Allowing  ior  -the- -false  gain  of  the  MCFs,  the  improvement  at 

_ _ - — - — -  ‘  ~  V 

higher  frequencies  (above  .75  Hz)  is  insignificant  and  a  simple  summation  pro- 

\  *»  1 

cess  is  considered  adequate.  The  improve 

band  is  directly  related  to  the  increase  in  the  surface  mode  noise  power  of  the 
noise  sample. 


rhent  obtained  in  the  .12  to  .75  Hz 


& 


C.  WAVENUMBER  ANALYSIS  OF  MCF  RESPONSES 

i 

The  wavenumber  response  of  the  MCFs,  MCF  W-34  and  MCF 
S-33,  were  computed  at'several  selected  frequencies  for  comparison  to  the 
wavenumber  spectra  of  the  noise  previously  discussed  in  Section  IV. 
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V-5 


services  group 


The  wavenumber  response  of  MCF  S-33,  at  the  lower  frequencies 
of  .  135  and  .  169  Hz  showing  the  k-space  area  being  significantly  rejected  is 
presented  in  Figure  V-3.  Response  at  the  signal  region  (K=0)  is  approximately 
.  0  db.  A  dashed  circle  in  the  spectra  represents  surface  mode  values.  Pri¬ 
marily  the  MCF  is  working  to  reject  the  surface  mode  energy  from  the  North¬ 
east  and  the  West-Southwest  with  some  emphasis  placed  on  rejecting  the  p-wave 
energy  coming  from  the  South.  There  are  also  indications  that  the  MCF  is 
attempting  to  reject  surface  mode  energy  from  the  direction  of  the  Gulf  of 
Mexico  to  the  Southeast  of  the  array. 

I 

At  two  higher  frequencies  of  .24  and  .30  Hz  (Figure  V-4)-the 
response  becomes  isotropic  to  surface  mode  velocities  with  some  emphasis  on 
the  energy  from  the  Southwest.  The  MCF  continues  to  reject  the  p-wave  energy 
from  the  South-Southwe3t.  The  velocity  of  16  km/second  from  the  Southwest 
is  strongly  Rejected  by  the  MCF  at  .  304  Hz. 

Comparable  wavenumber  responses  for  the  MCF  W-34  is  shown 
in  Figure  V-5.  At  .  122  Hz  the  MCF  is  strongly  rejecting  the  surface  mode 
energy  from  the  Northwest  and  Southeast  to  obtain  the  23  db  signal-to-noise 
improvement  of  Figure  V-l.  As  the  frequency  is  increased  to  ,244  Hz,  the 
MCF  begins  to  significantly  reject  the  p-wave  energy  from  the  Northwest  in 
addition  to  isotropic  rejection  of  surface  mode  energy.  At  higher  frequencies 
the  MCF  rejection  of  p-wave  energy  also  becomes  isotropic. 

The  significant  difference  in  the  'winter'  and  'summe^'  MCF  is 
the  rejection  of  surface  mode  energy  from  the  Northwest&t  low  frequencies 
(.  12  to  .24  Hz)  and  a  more  isotropic  p-wave  response  at  higher  frequencies 
(.  3  to  .  4  Hz). 
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D.  CONCLUSIONS 

The  performance  of  multichannel  filtering  for  the  extended 

short-period  array  is  highly  dependent  on  the  level  of  the  noise.  Under  noisy 

'j  I 

conditions,  apparently  caused  by  storm  activity  near  coastlines,  MCFs  can 
provide  significant  improvement  over  beemsteer  processing  at  frequencies 
below  ,75  Hz.  Above  this  frequency,  however,  very  little  improvement  is 

i  e 

obtained  under  either  condition.  Since  the  frequencies  of  prime  signal  interest 

In 

lie  in  the  higher  frequency  band  and  the  generating  sources,  below  .75  Hz  are 
not  stationary,  the  designed  MCFs  would  not  be  useful  for  on-line  application. 
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Several  contributors  to  the  ambient  noise  field  could  be  defined  in  both  of  the 
periods  analyzed  and  appeared  to  be  both  time-and  space  stationary.  Other  contri¬ 
butors,  due  primarily  to  atmospheric  storm  or  pressure  activity,  were  not  station¬ 
ary  and  could  be  related  to  it,s  gene-ra-ting-  source.  _  _  . 

Multi-channel  filtering  of  the  extended  array  could  not  provide  any  significant 
improvement  over  simple  beamsteering  at  frequencies  of  interest  (greater  than  .8 
Hz).  At  lower  frequencies,  the  amount  of  improvement  obtained  by  MCF  processing 
relative  to  beamsteering  was  directly  related  to  the  increase  in  noise  level  due  to  > 
•the  non- stationary  contributor  of  the  noise  field.  , 
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